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We present a simultaneous measurement of the B° lifetime t b o and B°B° oscillation frequency 
Anid- We use a sample of about 50 000 partially reconstructed B° — > D* + £~i>e decays identified 
with the BABAR detector at the PEP-II e + e~ storage ring at SLAC. The flavor of the other B 
meson in the event is determined from the charge of another high-momentum lepton. The results 
are 

t b o = (1.504 ± 0.013 (stat) to^l (syst)) ps, 
Am d = (0.511 ± 0.007 (stat) to'.ooe ( s y st )) P s_1 - 

PACS numbers: 13.25.Hw, 12.15.Hh, 14.40.Nd, 11.30.Er 
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I. INTRODUCTION 

The time evolution of B° mesons is governed by the 
overall decay rate T(B°) = 1/t b o and by the mass differ- 
ence Anid of the two mass eigenstates. A precise determi- 
nation of T(B°) reduces the systematic error on the pa- 
rameter \V c b\ of the Cabibbo-Kobayashi-Maskawa quark 
mixing matrix [1]. The parameter iVtdV^I enters the box 
diagram that is responsible for B°B° oscillations and can 
be determined from a measurement of Amd, although 
with sizable theoretical uncertainties. 

We present a measurement of t b o and Arrid using 
B a — > D* + l~i>£ decays [2] selected from a sample of 
about 88 million BB events recorded by the BABAR de- 
tector at the PEP-II asymmetric-energy e+e~ storage 
ring, operated at or near the T(4S) resonance. BB pairs 
from the T(4S) decay move along the beam axis with a 
nominal Lorentz boost (fl-f) — 0.56, so that the vertices 
from the two B decay points are separated on average by 
about 260 ^m. The B°B° system is produced in a coher- 
ent P-wave state, so that flavor oscillation is measurable 
only relative to the decay of the first B meson. Mixed 
(unmixed) events are selected by the observation of two 
equal (opposite) flavor B meson decays. The probabili- 
ties of observing mixed («S _ ) or unmixed (S + ) events as a 
function of the proper time difference At between decays 
are 

e -\At\/r B0 

S ± (At) = (l±X>cos(Am d At)), (1) 

4r B o 

where the dilution factor V is related to the fraction w of 
events with wrong flavor assignment by the relation V = 
1 — 2w and At is computed from the distance between 
the two vertices projected along the beam direction. 



II. THE BABAR DETECTOR AND DATASET 

We have analyzed a data sample of 81 fb _1 collected by 
BABAR on the T(45) resonance, a sample of 9.6 fb _1 col- 
lected 40 MeV below the resonance to study the contin- 
uum background, and a sample of simulated BB events 
corresponding to about three times the size of the data 
sample. The simulated events are processed through the 
same analysis chain as the real data. BABAR is a multi- 
purpose detector, described in detail in Ref. [3]. The mo- 
mentum of charged particles is measured by the tracking 
system, which consists of a silicon vertex tracker (SVT) 
and a drift chamber (DCH) in a 1.5-T magnetic field. 
The positions of points along the trajectories of charged 
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tracks measured with the SVT are used for vertex recon- 
struction and for measuring the momentum of charged 
particles, including those particles with low transverse 
momentum that do not reach the DCH due to bend- 
ing in the magnetic field. The energy loss in the SVT 
is used to discriminate low-momentum pions from elec- 
trons. Higher-energy electrons are identified from the 
ratio of the energy of their associated shower in the elec- 
tromagnetic calorimeter (EMC) to their momentum, the 
transverse profile of the shower, the energy loss in the 
DCH, and the information from the Chcrenkov detector 
(DIRC). The electron identification efficiency is about 
90%, and the hadron misidentification probability is less 
than 1%. Muons arc identified on the basis of the en- 
ergy deposited in the EMC and the penetration in the 
instrumented flux return (IFR) of the superconducting 
coil, which contains resistive plate chambers interspersed 
with iron. Muon candidates compatible with the kaon 
hypothesis in the DIRC are rejected. The muon identifi- 
cation efficiency is about 60%, and the hadron misiden- 
tification rate is about 2%. 



III. ANALYSIS METHOD 

A. Selection of B° -f D*+£~ v t decays 

We select events that have more than four charged 
tracks. We reduce the contamination from light-quark 
production in continuum events by requiring the normal- 
ized Fox- Wolfram second moment [4] to be less than 0.5. 
We select B° — ► D* + l~i>t events with partial reconstruc- 
tion of the decay D* + — ► n+D°, using only the charged 
lepton from the B° decay and the soft pion (7r+ ) from 
the D* + decay. The D° decay is not reconstructed, re- 
sulting in high selection efficiency. BABAR has already 
published two measurements of t b o [5, 6] and a measure- 
ment of sin(2/3 -(- 7) [7] based on partial reconstruction 
of B decays. This technique was originally applied to 
B° -> D*+£-D e decays by ARGUS [8], and then used by 
CLEO [9], DELPHI [10], and OPAL [11]. 

To suppress leptons from several background sources, 
we use only high-momentum leptons, in the range 1.3 < 
pe < 2.4 GeV/c [12]. The 7r+ candidates have momenta 
{p n +) between 60 and 200 MeV/c. Due to the limited 
phase space available in the D* + decay, the 7r+ is emitted 
within an approximately one-radian half-opening-angle 
cone centered about the D* + flight direction. We ap- 
proximate the direction of the D* + to be that of the 7r+ 
and estimate the energy E D *+ of the D* + as a function 
of the energy of the 7r+ using a third order polynomial, 
with parameters taken from the simulation. We define 
the square of the missing neutrino mass as 

M 2 V = (j§-E D , + -E e -y -(p D . + +p e -) 2 , (2) 

where we neglect the momentum of the B° in the Y(4S) 
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FIG. 1: Ml distribution for right-charge (top) and wrong- 
charge (bottom) events. The points correspond to on- 
rcsonance data. The distributions of continuum events (dark 
histogram), obtained from luminosity- rescaled off- resonance 
events, and BB combinatorial background events (hatched 
area), obtained from the simulation, are overlaid. Monte 
Carlo events are normalized to the difference between on-peak 
and rescaled off peak data in the region Ml < —4.5 GeV 2 /c 4 . 



frame (on average, 0.34 GeV/c), and identify the B° en- 
ergy with the beam energy in the e + e~ center-of- 
mass frame. Eg- and p e are the energy and momentum 
vector of the lepton and is the estimated momen- 
tum vector of the D* + . The distribution of Ml peaks 
at zero for signal events, while it is spread over a wide 
range for background events (see Fig. 1). 

We determine the B° decay point from a vertex fit of 
the l~ and 7r+ tracks, constrained to the beam-spot posi- 
tion in the plane perpendicular to the beam axis (the x-y 
plane). The beam-spot position and size are determined 
on a run- by-run basis using two-prong events [3] . Its size 
in the horizontal (x) direction is on average 120 /im. Al- 
though the beam-spot size in the vertical (y) direction is 
only 5.6 /im, we use a constraint of 50 /xm in the vertex 
fit to account for the flight distance of the B° in the x — y 
plane. We reject events for which the x 2 probability of 
the vertex fit, Vv, is less than 0.1%. 



We then apply a selection criterion to a combined sig- 
nal likelihood, X, calculated from p^-, p +, and Vv, 
which results in a signal-to-background ratio of about 
one in the signal region defined as Ml > —2.5 GeV 2 /c 4 . 
We reject events for which X is lower than 0.4 (see Fig. 2). 
Figure 1 shows the distribution oi Ml after this selec- 
tion. The distributions in the top part of the figure are 
obtained from events in which the I and the n s have op- 
posite charges ("right-charge"), and the distributions in 
the bottom are from events in which the t and the tt s 
have equal charges ( "wrong-charge" ) . 

The points in Fig. 1 correspond to on-resonance data. 
The dark histograms correspond to off-resonance data, 
scaled by the ratio of on-rcsonancc to off-resonance in- 
tegrated luminosity. The hatched histograms correspond 
to BB combinatorial background from simulation. To 
normalize the BB combinatorial background, we scale 
the BB Monte Carlo histogram so that, when added to 
the luminosity-scaled off-resonance histogram, the sum 
matches the on-resonance data in the region Ml < —4.5 
GeV 2 /c 4 . The right-charge plot is shown for illustra- 
tion only. We use the wrong-charge samples as a cross- 
check to verify that the BB combinatorial background 
shape is described by the simulation. For this purpose, 
we compare the number of wrong-charge events in the sig- 
nal region predicted from the sum of off-resonance and 
BB Monte Carlo, normalized as above, to the number 
of wrong-charge on-resonance data events. This ratio is 
0.996 ± 0.002, consistent with unity. For the rest of the 
analysis we consider only right-charge events. 



B. Tag Vertex and B Flavor Tagging 

To measure Amj we need to know the flavor of both B 
mesons at their time of decay and their proper decay time 
difference At. The flavor of the partially reconstructed 
B is determined from the charge of the high- momentum 
lepton. In order to identify the flavor of the other ("tag") 
B meson, we restrict the analysis to events in which an- 
other charged lepton (the "tagging lepton") is found. To 
reduce contamination from fake leptons and leptons orig- 
inating from charm decays, we require that the momen- 
tum of this second lepton exceed 1.0 GeV/c for electrons, 
and 1.1 GeV/c for muons. 

The decay point of the tag B is determined with the 

high- momentum lepton and a beam-spot constraint; the 

— o 

procedure is the same as that used to determine the B 
vertex. We compute At from the projected distance be- 
tween the two vertices along the beam direction (z-axis), 
Az = Zdecay — z ta g, with the approximation that the 
B a and the B° are at rest in the T(4S') rest frame (the 
boost approximation): At = Az/cf3-f, where the boost 
factor (3j is determined from the measured beam ener- 
gies. To remove badly reconstructed vertices we reject 
all events with either \Az\ > 3 mm or cr(Az) > 0.5 mm, 
where a(Az) is the uncertainty on Az, computed for each 
event. The simulation shows that the difference between 
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FIG. 2: Distribution of the combined signal likelihood X for 
events in the signal M v 2 region. Events for which X < 0.4 
are rejected. 



the true and measured At can be fitted with the sum of 
two Gaussians. The rms of the narrow Gaussian, which 
describes 70% of the events, is 0.64 ps; the rms of the 
wide one is about 1.7 ps. 

We then select the best right-charge candidate in each 
event according to the following procedure: if there is 
more than one, we choose that with Ml > —2.5 Gev 2 /c 4 . 
If two or more candidates are still left, but they have dif- 
ferent leptons, we select the one with the largest value 
of x- I n a small fraction of events we select two or 
more candidates sharing the same lepton combined with 
different soft pions. We keep the candidate with the 
largest x, unless one of the 7r+ is consistent with coming 
from the decay of a D* + from the other B, in which 
case we remove the event. For this purpose, we define 
the square of the missing neutrino mass in the tag-side, 
Ml tag _ B , by means of Eq. 2, where we replace the four- 
momentum of the lepton from the B° — ^ D* + £~vi decay 
with that of the tag lepton. This variable peaks at zero 
for soft pions originating from the tag-B decay. We re- 
quire Ml tag _ B < —3 GeV 2 /c 4 . Finally we reject the 
events in which the signal lepton can be combined to 
a wrong-charge pion to produce an otherwise successful 
candidate, if the pion is consistent with coming from a 



D* + from the tag-B according to the criterion just de- 
scribed. About 20% of the signal events are removed by 
this requirement. 

For background studies, we select events in the region 
Ml < —2.5 GeV 2 /c 4 if there is no candidate in the signal 
region. We find about 49000 signal events over a back- 
ground of about 28000 events in the data sample in the 
region Ml > -2.5 GeV 2 /c 4 . 

C. Sample Composition 

Our data sample consists of the following event types, 
categorized according to their origin and to whether or 
not they peak in the Ml distribution. We consider sig- 
nal to be any combination of a lepton and a charged 
D* produced in the decay of a single B° meson. Sig- 
nal consists of mainly B° — > D* + £~V£ decays, with minor 
contributions from B° -> D* + Tfrv t , B° -> D*+t_V t , 
B° D*+D-, and B° D^+DX with r, D s ~ , or D de- 
caying to an £~ , and from B° — ► D* + h, with the hadron 
h misidentified as a muon. Peaking B~ background is 
mainly due to the processes B~ — > D* + ir~ £~Vi, and 
B~ — > D* + tt~X with the 7r~ misidentified as a muon. 
Other minor contributions to the peaking sample are 
due to decays B D*txv t t (t -> tX), B D*ttDX 
(D — > £Y), where the D* and the it come from the decay 
of an orbitally excited D meson (D**). Non-peaking con- 
tributions are due to random combinations of a charged 
lepton candidate and a low-momentum pion candidate, 
produced either in BB events (BB combinatorial) or in 
e + e _ — > qq interactions with q = u, d, s, or c (contin- 
uum). We compute the sample composition separately 
for mixed and unmixed events by fitting the correspond- 
ing Ml distribution to the sum of four components: con- 
tinuum, BB combinatorial background, B° — * D* + i~h*i 
decays, and B — > D*it£~vi decays. Due to one or more 
additional pions in the final state, the B — > D*nl~vg 
events have a different M v 2 spectrum from that of the 
process B° — > D* + £~i/£. We measure the continuum con- 
tribution from the off-resonance sample, scaled to the 
luminosity of the on-resonance sample. We determine 
the Ml distributions for the other event types from the 
simulation, and determine their relative abundance in 
the selected sample from a fit to the Ml distribution 
for the data. Assuming isospin conservation, we assign 
two thirds of B — > D*ir£~i7g decays to peaking B~ back- 
ground and the rest to £>° — > D*tt£~Vi , which we add to 
the signal. We vary this fraction in the study of system- 
atic uncertainties. We assume 50% uncertainty on the 
isospin-conservation hypothesis. 

A possible distortion in the A^ 2 distribution comes 
from the decay chain B — > D(X)£~~p£, D — > Yir + , where 
the state Y is so heavy that the charged pion is emitted at 
low momentum, behaving like a 7r+ . This possibility has 
been extensively studied by the CLEO collaboration [13], 
where the three D + decay modes most likely to cause this 
distortion have been identified: K W7r + , K*~ p + ir + , and 
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K p 7r + . If wc remove these events from the simulated 
sample, and we repeat the fit, the number of fitted signal 
events increases by 0.4%. We assume therefore ±0.4% 
systematic error on the fraction of signal events in the 
sample due to this uncertainty. 

Figure 3 shows the Ai v 2 fit results for unmixed (up- 
per) and mixed (lower) events. We use the results of this 
study to determine the fraction of continuum (f^ q ), BB 
combinatorial (f^-), and peaking B~ {fg- ) background 
as a function of A4 2 ,, separately for mixed (/ _ ) and un- 
mixed (/ + ) events. We parameterize these fractions with 
polynomial functions of Mj 2 as shown in Fig. 4. 



D. t b o and Amd Determination 

We fit data and Monte Carlo events with a binned 
maximum-likelihood method. We divide the events into 
one hundred At bins, spanning the range — 18 ps < At < 
18 ps, and twenty cta* bins between and 3 ps. We as- 
sign to all events in each bin the values of At and ca* 
corresponding to the center of the bin. We fit simultane- 
ously the mixed and unmixed events. We maximize the 
likelihood 



n 




a 



s„- x 



(3) 



where 
mixed 



the indices k and j denote the un- 
and mixed selected events. The functions 
J-^ (At, aAt, -MUtbo , Amd) describe the normalized At 
distribution as the sum of the decay probabilities for 
signal and background events: 



J 7± (At 1 £TA 1 ,^|r B o,Am (i ) 

ft q (Ml)-J^ g (At,a At ) 



qq\ 
7± 
BB 1 



(4) 



•F^(At,<TAi) 



+S B -f±4Ml)-F±_(At,a A t) 
+ [1 - S B -fa(Ml) - fi B {Ml) - f±(Ml)] 
■Fgo (A*, VAt\TB» , Am d ), 

where the functions Tf~ represent the probability den- 
sity functions (PDF) for signal (i = B°), peaking _B~ 
(i = B~), BB combinatorial (i = BB), and continuum 
(i = qq) events, modified to account for the finite resolu- 
tion of the detector, and the superscript +(— ) applies to 
unmixed (mixed) events. The resolution function is ex- 
pressed as the sum of three Gaussian functions, described 
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FIG. 3: Fit to the M v distribution for the unmixed events 
(a) and mixed events (b). includes B° — > D* + l~vi, 

B° -» D* + ir°r_v e , B° -> ZT+r j? T (r -> IX), B^_-+ D*+ D~ 
[D s -> IX), B° -» D*+DX (D -> £Y), and B° -> D*+h 
with the hadron h misidentified as a muon. includes 
B~ -> D* + TvTu e and B~ -> D* + ir~X with the tt" 
misidentified as a muon. "Other" includes B — » D*m> T T 
(t -> IX), B -> D'tvDX (D -> IY). 
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as "narrow" , "wide" , and "outlier" : 



K(8At,(r A t) 



\/27rS n aAt 

fw 

\/2nS w <TAt 

f (i5At-o ) 2 
Jo _ 7^2 



(sAt— o w y 



where SAt is the difference between the measured and 
true values of At, o n and o w are offsets, and the factors 
S n and S w account for possible misestimation of a At- 
The outlier term, described by a Gaussian function of 
fixed width S and offset o a , is introduced to describe 
events with badly measured At, and accounts for less 
than 1% of the events. 

To account for the ±50% uncertainty on the isospin as- 
sumption (see Sec. IIIC), the functions fg_ and are 
multiplied in the PDF for the peaking B~ background 
by a common scale factor Sg - . This parameter is al- 
lowed to vary in the fit, constrained to unity with vari- 
ance erg- = 0.5 by means of the Gaussian term 



(5 F 



a 



We constrain the expected fraction P cxp 
events to the observed one 



of mixed 



Pahs = 



•A/mix 



Mr 



■M u: 



by means of the binomial factor 



M! 



(1 = 
Xd Af ■ W 



T exp 



'(1 -Pexp) 



For a sample of signal events with dilution T>, the ex- 
pected fraction reads 



P e xp(Amrf, T B o , V) = Xd ■ V 



1 - V 



where, neglecting the decay-rate difference Ar^ between 
the two mass eigenstates, the integrated mixing rate Xd 
is related to the product x = Am c i ■ Tgo by the relation 



FIG. 4: Fraction of continuum events, peaking B~ , and BB 
combinatorial background in the unmixed (a) and mixed (b) 
lepton-tagged samples. The continuous lines overlaid repre- 
sent the analytic functions (/^, and fgg) used to param- 
eterize the distributions. The fraction of continuum events is 
parameterized only in the region M„ > —4.5 GeV 2 /c 4 . For 
J\A„ < —4.5 GeV 2 /c 4 , where just continuum and combinato- 
rial backgrounds are present, we assume that ftj_ = 0, and 



we compute fq 



f 



Xd 



2(1 + x 2 ) 



We divide signal events according to the origin of the 
tag lepton into primary (VI), cascade (C£), and decay- 
side (D£) lepton tags. A primary lepton tag is produced 
in the direct decay B° — > l + vgX. These events are de- 
scribed by Eq. 1, with V close to 1 (a small deviation 
from unity is expected due to hadron misidentification, 
lcptons from J/tjj, etc.). We expect small values of o„ 
and o w for primary tags, because the lepton originates 
from the B° decay point. 

Cascade lepton tags, produced in the process B° — > 
DX, D — > £Y, are suppressed by the requirement on the 
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lepton momentum but still exist at a level of 9%, which 
we determine by varying their relative abundance (fee) 
as an additional parameter in the Am^ and t^o fit 
on data. The cascade lepton production point is dis- 
placed from the B° decay point due to the finite lifetime 
of charm mesons and the e + e~ energy asymmetry. This 
results in a significant negative value of the offsets for this 
category. Compared with the primary lepton tag, the 
cascade lepton is more likely to have the opposite charge 
correlation with the B° flavor. The same charge corre- 
lation is obtained when the charm meson is produced 
from the hadronization of the virtual W from B decay, 
which can result in the production of two opposite-flavor 
charm mesons. We account for these facts by applying 
Eq. 1 to the cascade tag events with negative dilution 

V ce = -(1 - 2flj- c ^ i+ ) = -0.65 ± 0.08, where we take 
from the PDG [14] the ratio: 



,6-> 

Ice 



B(b 



B(b -> c -> £+) + B(b^c^ i-) 



0.17 ±0.04. 



The contribution to the dilution from other sources as- 
sociated with the tt+£~ candidate, such as fake hadrons, 
is negligible. 

Decay-side tags are produced by the semileptonic de- 
cay of the unreconstructed D°. Therefore they do not 
carry any information about Tgo or Arra^. The PDF for 
both mixed and unmixed contributions is a purely expo- 
nential function, with an effective lifetime td c represent- 
ing the displacement of the lepton production point from 
the B° decay point due to the finite lifetime of the D°. 
We determine the fraction of these events by fitting the 
cos0 v + t - distribution (see plots in Figs. 5 (a) and 6 (a)), 
where is the angle between the soft pion and the 

tag lepton in the e + e~ rest frame. We fit the data with 
the sum of the histograms for signal events, BB com- 
binatorial background, and peaking B~ background ob- 
tained from the simulation, and continuum background 
obtained from the off-resonance events. We fix the frac- 
tion of signal events, peaking B~ background, BB com- 
binatorial background and continuum background in the 
fit and we allow to vary the relative amount of decay-side 
tags and tag-side tags. 

Using the results of the cos + e - ^ wc parameterize 
the probability for each event to have a decay-side tag as 
a third-order polynomial function of cos O^+e- ( see plots 
in Figs. 5 (b) and 6 (b)). 

The signal PDF for both mixed and unmixed events 
consists of the sum of PDFs for primary, cascade, and 
decay-side tags, each convoluted with its own resolution 
function. The parameters S n , S w , S a , f w , and f a are 
common to the three terms, but each tag type has differ- 
ent offsets (o„, o w ). All the parameters of the resolution 
functions, the dilution of the primary tags, the fraction 
of cascade tags, and the effective lifetime of the decay- 
side tags are free parameters in the fit. We fix the other 
parameters (dilution of cascade tags, fraction of decay- 
side tags) to the values obtained as described above, and 
then vary them within their uncertainties to assess the 
corresponding systematic error. 
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FIG. 5: Distribution of cos 9 +,_ for unmixed events. In fig- 
ure (a), the points with error bars represent the data, and the 
histograms show the various sample components determined 
by the fit. Figure (b) shows the ratio between the fraction of 
tags from D° decays over the total number of tags in signal 
events, as obtained from the simulation rescaled to the result 
of the fit. This distribution is parametrized by a third order 
polynomial represented by the continuous line overlaid. The 
dotted line represents the corresponding distribution obtained 
without rescaling the simulation to the fit result. 
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FIG. 6: Same as Fig. 5 for mixed events. 



We adopt a similar PDF for peaking B~ background, 
with separate primary, cascade, and decay-side terms. 
Because B~ mesons do not oscillate, we use a pure ex- 
ponential PDF for the primary and cascade tags with 
lifetime t#- = 1.671 ps [15]. We force the parameters of 
the resolution function to equal those for the correspond- 
ing signal term. 

We describe continuum events with an exponential 
function convoluted with a three-Gaussian resolution 



function. The mixed and unmixed terms have a common 
effective lifetime Ti q . All the parameters of the continuum 
resolution function are set equal to those of the signal, 
except for the offsets, which are free in the fit. 

The PDF for combinatorial BB background accounts 
for oscillating and non-oscillating subsamples. It has the 
same functional form as the PDF for peaking events, 
but with independent parameters for the oscillation fre- 
quency, the lifetimes and the fractions of B~ background, 
primary, cascade, and decay-side tag events. The param- 
eters S w , f w and o are set to the same values as those 
in the signal PDF. 



IV. RESULTS 

We first apply the measurement procedure on several 
Monte Carlo samples. We validate each term of the 
PDF by first fitting signal events, for primary, cascade 
and decay-side tags separately, and then adding them 
together. We then add peaking B~ background, and fi- 
nally add the BB combinatorial background. We observe 
the following features: 

• The event selection introduces no bias on t b o and 
a bias of (-0.0029 ± 0.0010) ps -1 on Am d . 

• The boost approximation introduces a bias on 
t b o (+0.0054 ps) and an additional bias on Am d 
(—0.0034 ps -1 ), determined by fitting the true Az 
distribution. These biases disappear however when 
we fit the smeared Az and allow for the experimen- 
tal resolution in the fit function. 

• After the introduction of B~ peaking background 
we observe a bias of (—0.0079 ± 0.0064) ps on r B a 
and (-0.0034 ±0.0028) ps" 1 on Am d . 

• Adding combinatorial BB events, we observe a bias 
of (+0.0063 ± 0.0070) ps on t b o and (-0.0074 ± 
0.0035) ps" 1 on Am d . 

• The isospin scale factor S B - = 0.91 + 0.10 is con- 
sistent with unity. 

Based on these observations, we correct the data re- 
sults by subtracting 0.0063 ps from t b o , and adding 
0.0074 ps -1 to Am d . We include the Monte Carlo statis- 
tical errors of ±0.0070 ps for t b o and ±0.0035 ps -1 for 
Amd as systematic uncertainties. 

We determine the parameters for continuum events di- 
rectly from the fit to on-resonance data, and we indepen- 
dently fit the off-resonance events to verify the consis- 
tency with the on-resonance continuum results. 

We finally perform the fit to the on-resonance data. 
Together with Am d and t b o, we allow to vary most of 
the parameters describing the peaking B~ , BB combina- 
torial, and continuum background events. The results of 
the fits to the Monte Carlo and data samples are shown 
in Tables I and II. 
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The fit results are t b o = (1.510 ±0.013 (stat)) ps, and 
Am d = (0.5035 ± 0.0068 (stat)) ps -1 . We correct these 
values for the biases measured in the Monte Carlo simu- 
lation, obtaining the results 

t b o = (1.504 ±0.013 (stat)) ps, 
Am d = (0.5109 ± 0.0068 (stat)) ps -1 . 

The statistical correlation between Am d and t b o is 0.7%. 
Amd has sizable correlations with S B - (50%) and with 
the fraction of cascade tags (24%). t b o is correlated with 
S n (—27%) and the offset of the wide Gaussian for the 
cascade tags (—31%). The complete set of fit parameters 
is reported in Tables I and II. 

Details on the systematic error are reported in Sec. V. 
Figures 7 and 8 show the comparison between the data 
and the fit function projected on At, for a sample of 
events enriched in signal by the cut Ml > —2.5 GeV 2 /c 4 ; 
Figs. 9 and 10 show the same comparison for events in 
the background region. 

Figures 11 and 12 show plots of the time-dependent 
asymmetry 

A(At) _ MmmixCAt) - A/"mix(Af) 
~ A/" un mix(At) + M mix (At) 

for events in the M. 2 signal region and events in the A4 2 , 
background region. For signal events, neglecting At res- 
olution, A(At) = V cos(Amd At) (see Eq. 1). 

The agreement between the fit function and the data 
distribution is good in both the signal and background 
regions. The asymmetry is quite significant for events 
in the background A4 2 , region because a large fraction of 
these events are due to combinatorial B°B° background. 



V. SYSTEMATIC UNCERTAINTIES 




At(ps) 



FIG. 7: Distribution of At for unmixed (a) and mixed (b) 
events in the signal M„ region. The points show the data, 
the curve is the projection of the fit result, and the shaded 
areas from bottom to top are the contributions from contin- 
uum background, BB combinatorial background, peaking B~ 
background with decay-side tag, peaking B~ background with 
primary tag, signal with decay-side tag, signal with cascade 
tag, and signal with primary tag. 



The systematic errors are summarized in Table III. We 
consider the following sources of systematic uncertainty: 

1. Sample composition: We calculate a total uncer- 
tainty of ±1.3% on the number of signal events. 
This uncertainty is the quadratic sum of the statis- 
tical error in the M„ 2 fit (±1.2%), the systematic 
uncertainty on the shape of BB combinatoric back- 
ground from the test on the "wrong-charge" sample 
(±0.2%) (see sec. Ill A), and the additional system- 
atic uncertainty due to low-momentum pions from 
D+ decays (±0.4%) (see sec. IIIC). 

2. Analysis bias (entry b): We use the statistical error 
on the bias observed in the fit on the Monte Carlo 
sample. 

3. Signal and background PDF description: Most of 
the parameters in the PDF are free in the fit and 
therefore do not contribute to the systematic error. 



We vary the parameters that are fixed in the fit by 
their uncertainty, repeat the fit, and use the corre- 
sponding variation in t b o and Am d as systematic 
errors. We take the uncertainty on t b - (entry c), 
and on D cl (entry d) from the PDG [14]. We find 
that four parameters used in the description of the 
combinatorial background, as determined by the fit 
on the Monte Carlo sample, are not in agreement 
with the Monte Carlo truth. They are the fraction 
of cascade tag-side leptons in the unmixed event 
sample, fc e K u G , the fraction of decay-side tags in 
the mixed B° and the B~ event samples, and 
a B— m> respectively, and an additional parameter 
used in the description of the shape of the proper 
time difference At of the decay-side tagged mixed 
sample, f-pe,2- Therefore we fix them to the Monte 
Carlo prediction. We vary the value of each of them 
by 20% to compute the systematic error from the 
comparison with the default result, and we sum the 
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TABLE I: Parameters used in the PDFs. The upper set of parameters refers to peaking events; the lower one refers to those 
parameters of the resolution function that are common to all the event types. The second column shows how the parameters 
are treated in the fit. The third (fourth) column gives the result of the fit on data (MC) for free parameters and the value 
employed for the parameters that are fixed or used as a constraint. The quoted error is the statistical uncertainty from the fit 
for free parameters and the range of variation used in the systematic error determination for the others. The last column shows 
the sample in which the parameter is used. VI, CI and VI refer to primary, cascade and decay-side lepton tags, respectively. 
The parameters o, S, and / correspond to offsets, scale factors, and fractions in the resolution function. 



Parameter 


Usage 


data 


ML 


sample 


T B (ps) 


free 


1.510±0.013 


1.554±0.007 


B° VI, CI 


Amd(ps _1 ) 


free 


0.503±0.007 


0.465±0.004 


B° VI, Ct 


T De (PS) 


free 


0.12±0.04 


0.21 ±0.02 


B°, B-, 










bb vi 


T B - (ps) 


fixed 


1.671±0.018 


1.65 


B~ , BB 


S B - 


constr. 


1.05±0.15 


0.91±0.10 






(1.0±0.5) 








fee 


free 


0.095±0.006 


0.066±0.004 


B° 


D-pe 


free 


1.00±0.02 


0.970±0.006 


B°,V£ 


D C e 


fixed 


-0.65±0.08 


-0.545 


B°,C£ 


Ovt,n 


free 


-0.019±0.011 


-0.012±0.007 


B°, B~ VI 


(= 0-pl,w) 










oee,n 


free 


-0.18±0.07 


-0.43±0.07 


b° ce 


oce,w 


free 


2.8±1.1 


-5.8±0.7 


B° CI 


o-ve,n 


free 


-0.12±0.03 


-0.14±0.02 


B°, B~, VI 


(= ox>e,w) 










S n 


free 


0.952±0.015 


1.007±0.006 


all 


So (ps) 


free 


12.8±5.6 


17.9±8.2 


B°, B~ 


fo 


free 


0.0013±0.0005 


0.0008±0.0003 


B°, fl- 


Sw 


free 


2.57±0.13 


2.63±0.15 


ail 


ftv 


free 


0.050±0.005 


0.035±0.005 


all 


Oo 


fixed 








all 



four uncertainties in quadrature (entry e). 

4. Detector alignment: We consider effects due to the 
detector z scale, determined by reconstructing pro- 
tons scattered from the beam pipe and compar- 
ing the measured beam pipe dimensions with the 
optical survey data [16]. The z scale indetermina- 
tion corresponds to an uncertainty of ±0.4% on At. 
We repeat the fit applying this scale correction to 
At, and use the variation with respect to the de- 
fault result as the systematic error (entry f). From 
the measurement of the beam energies, the T(4S) 
Lorentz boost factor is determined with an uncer- 
tainty which translates into a ±0.1% indetermina- 
tion on At. Again we repeat the fit and assume 
as systematic error the variation of the result (en- 
try g). We then consider the effect of varying the 
beam-spot position by ±40 fim in the y direction 
(entry h). We compute the uncertainty due to SVT 
time-dependent misalignment by comparing results 
obtained with different sets of alignment constants 
(entry i). 

5. Decay-side tags: We vary the parameters describ- 
ing the fraction of decay-side tags by their statis- 
tical errors, repeat the fit, and take the variation 



with respect to the default result as the systematic 
error (entry j). 

6. Binned fitting: We vary the number of bins in At 
from 100 to 250 and in aAt from 20 to 50, and 
we repeat the fit. We take the systematic error to 
be the variation with respect to the default result 
(entry k). 

7. Outlier description: We vary the value of the offset 
of the outlier Gaussian from —5 ps to +5 ps. As a 
cross-check, we use a PDF that is uniform in At for 
the description of the outliers. We take the maxi- 
mum variation with respect to the default result as 
the systematic error (entry 1). 

8. Fit range: We vary the At fit range from ±18 ps to 
± 10 ps and the a At maximum value from 1.8 ps to 
4.2 ps. Again we assume the maximum variation 
between the various results and the default one as 
the systematic error (entry m). 

9. Cascade lepton tag-side parameterization: For the 
resolution model, we use a Gaussian distribution 
convolved with a one-sided exponential to describe 
the core part of the resolution function (GExp) in- 
stead of the Gaussian resolution with a non-zero 
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TABLE II: Parameters used in the background PDF. The upper set of parameters refers to BB combinatorial events, the 
central one refers to continuum parameters, and the lower set refers to those parameters of the resolution function that are 
common to all event types. The symbols a correspond to the fractions of decay-side tags in the different samples. The last line 
shows the statistical correlation between t b o and Arrid. 



Parameter 


Usage 


data 


M.C. 


Sample 


T B B KG (PS) 

Am B d KG (ps- 1 ) 
r B - BKG (ps) 

rBKG 

Jce,u 

rBKG 
JCi,m 

rBKG 
JB-,u 

rBKG 
JB-,m 
BKG 
a BO.u 
BKG 
a BO.m 

&B— u 
BKG 

a B-,m 

fae,2 

j-\BKG 

nBKG 
u Cl 
BKG 

°ve 

BKG 
Oct 

qBKG 
>~>n 

S BKG (ps) 

rBKG 
Jo 


free 
free 
fixed 
fixed 
free 
free 
free 
free 
fixed 
free 
fixed 
fixed 
free 
fixed 
free 
free 
free 
free 
free 
free 


1.22±0.07 
0.37±0.06 
1.671±0.018 
0.030±0.006 
0.041±0.022 
0.62±0.08 
0.15±0.10 
0.11±0.19 
0.065±0.013 
0.21±0.10 
0.36±0.07 
0.60±.12 
0.989±0.013 
-0.65±0.08 
-0.006±0.016 

— l.OztU.D 

-0.02±0.04 
0.961±0.015 
10.4±3.6 
0.0021±0.0009 


1.37±0.07 
0.42±0.04 
1.65 
0.030 
0.069±0.021 
0.52±0.02 
0.11±0.04 
0.25±0.03 

0.065 
0.20±0.02 

0.36 
0.60±.12 
0.964±0.006 

-0.545 
-0.02±0.03 

H £-1-0 9 

-0.05±0.03 
0.961±0.021 
14.7±6.1 
0.0008±0.0003 


BB 

BB 
B~, BB 
BB (B° only) 
BB (B° only) 

BB 

BB 
BB (B° only) 
BB (B° only) 
BB (B~ only) 
BB (B~ only) 

BB, VI 
BB (W only) 
BB {B° only) 
BB, VI 

BB, VI 
all 
BB 
BB 


Tlq (PS) 
Olq.n 

(= Olq }W ) 


free 
free 


0.27±0.05 
0.007±0.032 




Continuum 
Continuum 


S w 
fw 
Oo 


free 
free 
fixed 


2.57±0.13 
0.050±0.005 



2.63±0.15 
0.035±0.005 



all 
all 
all 


p( T B a ,Am d ) 




0.007 


-0.127 





offset. We quote as systematic error the difference 
between the results obtained with the two different 
approaches (entry n). 



VI. CONSISTENCY CHECKS 

We rely on the assumption that the parameters of 
the background PDF do not depend on Ml- We verify 
this assumption for the continuum background with the 
fit to the off-resonance events. To check this assump- 
tion for the BB combinatorial PDF, we perform sev- 
eral cross checks on the data and the Monte Carlo. We 
compare the simulated combinatorial BB At distribution 
in several independent regions of Aif, with Kolmogorov- 
Smirnov tests and always obtain a reasonable probability 
for agreement. We fit the At distribution of combinato- 
rial background BB events separately in the signal and 
background ]\A 2 V region and compare the parameters of 
the PDF. We fit the signal plus BB background Monte 
Carlo events in the signal region only, fixing all the pa- 
rameters of the BB background to the values obtained in 



a fit to the background region, and do not see any signif- 
icant deviation from the results of the full fit. Finally, we 
repeat the fit on both the data and the Monte Carlo us- 
ing different M 2 V ranges for the background region. Once 
again, we do not observe any significant difference in t b o 
and Arrid relative to the default result. 

We repeat the analysis with a more stringent require- 
ment on the combined signal likelihood (a minimum X of 
0.5 rather than 0.4). No significant change in the result 
is observed. 

We validate the fit procedure with a parameterized 
Monte Carlo simulation. We simulate several experi- 
ments from the fitted PDF of both the Monte Carlo and 
the data, with parameters fixed to the values obtained 
from the corresponding fit. Each experiment is produced 
with the same number of events as the original sample. 
For each experiment we produce seven data sets, corre- 
sponding to B° with primary, cascade, and decay-side 
lepton tags, peaking B~ background with tag-side and 
decay-side lepton tags, BB combinatorial background, 
and continuum background. We fit every experiment 
with the same procedure as the corresponding original 
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FIG. 8: Same as Fig. 7 with logarithmic scale. 



TABLE III: Systematic uncertainties. See text for details. 



Source 



(a) Sample Composition 

(b) Analysis bias 

( c ) T B- 

(d) Vci 

(e) Combinatorial BKG 

(f) z scale 

(g) PEP-II boost 

(h) Beam-spot position 

(i) Alignment 

(j) Decay-side tags 
(k) Binning 
(1) Outlier parameters 
(m) At and a At cut 
(n) GExp model 



Total 



Variation 



±1.3% 

1.671±0.018 
0.65±0.08 



8r B o (ps) SArrid (ps 



±0.0003 
±0.0070 
+"0.0014 
+0.0003 
±0.0007 
±0.0070 
±0.0020 
±0.0050 

+0.0132 
-0.0038 

±0.0013 
+0.0021 
±0.0028 
±0.0033 
-0.0016 



+0.0182 

-0.0131 



+0.0002 
+0.0035 
+0.0008 
+0.0003 
+0.0002 
+0.0020 
+0.0003 

+0.0010 

-0.0038 
+0.0033 

±0.0006 
±0.0012 
+0.0033 
±0.0011 



+0.0068 
-0.0064 



FIG. 9: Distribution of At for unmixed (a) and mixed (b) 
events in the background M 2 V region. The points show the 
data, the curve is the projection of the fit result, and the 
shaded areas are the contributions from continuum and BB 
combinatorial background. 



We do not find any significant statistical anomaly in the 
fit behaviour. 



TABLE IV: Results extracted from parameterized Monte 
Carlo experiments generated with parameters fixed to the val- 
ues obtained from the fit to data (second column) and to the 
full Monte Carlo simulation (third column). For both r B o and 
Arrid the average and the rms of the distribution of the pull 
with respect to the generated value are reported. 



Parameter 


data 


Monte Carlo 


number of experiments 
Pull r B a average 
Pull t b o rms 
Pull Arnj average 
Pull Arrid rms 


54 

0.38 ±0.19 

1.13 ±0.19 
-0.33 ±0.17 

1.14 ±0.17 


124 
0.38 ±0.12 
1.25 ±0.13 
0.08 ±0.11 
1.09 ±0.08 



sample, and finally we compare the fitted parameters 
with the generated values. The result of this study is 
summarized in Table IV where we report the average and 
the root-mean-square deviation (rms) of the distribution 
of the difference between the fitted and the generated pa- 
rameter value divided by the fit statistical error (pull). 



We rely on the assumption that the decay-rate differ- 
ence ATd between the two mass eigenstates can be ne- 
glected in the analysis. We check this assumption with 
a parameterized Monte Carlo simulation in which events 
are simulated with zero mistag probability and perfect At 



17 



10 ; 



U5 

d 



•S 10' 



10 





1 




10 3 


& 




vo 




CO 




o 






10 2 






•e 




a 






10 




1 




At(ps) 



FIG. 10: Same as Fig. 9 with logarithmic scale. 



resolution. We produce two sets of one hundred Monte 
Carlo experiments. In the first set, AI^ = 0; in the sec- 
ond, AJ-^ = 0.01. We fit every experiment with the same 
procedure neglecting Al^ and we do not find any signif- 
icant difference in the values of t^o and Amj in the two 
different sets. 

We investigate a possible analysis bias due to the fi- 
nite -r_and D s lifetimes in B° -> D*+t~v t (t~ -> i~ X) 
and B° -> D*+D- (D~ -► i~X) decays. We fit the 
Monte Carlo signal sample with no mistag and realistic 
At resolution after removing these decays and we do not 
find any significant variation with respect to the result 
obtained with the full signal sample. 



VII. CONCLUSION 

We have performed a measurement of Amd and t b o 
with a sample of about 50 000 partially reconstructed, 
lepton-tagged B° — > D* + £~Di decays. We obtain the fol- 
lowing results: 



T B o 

Am d 



(1.504 ±0.013 (stat) ±olox3 ( s y st )) P s ' 
(0.511 ±0.007 (stat) (syst)) ps" 



>4-5 
— 



< 



0.5 



-0.5 




-10 



10 



At(ps) 



FIG. 11: Asymmetry between unmixed and mixed events as a 
function of At, for events in the signal M„ region. Points with 
error bars represent the data, and the curve is a projection of 
the fit result. 



The Tgo value is consistent with the published measure- 
ment performed by BABAR using B Q — ► D* + £~P( par- 
tially reconstructed decays [5]. Our results are also con- 
sistent with published measurements of Tgo and Am<j 
performed by BABAR with different data sets [6, 17-20], 
and with the world averages computed by the Heavy 
Flavor Averaging Group for the PDG 2005 web up- 
date: t b o = (1.532 ± 0.009) ps, and Am d = (0.505 ± 
0.005) ps" 1 . 
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projection of the fit result. 
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